Coexistence of superconductivity and a spin density wave in pnictides: Gap symmetry and nodal 
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We investigate the effect of a spin-density wave (SDW) on s± superconductivity in Fe-based superconductors. 
We show that, contrary to the common wisdom, no nodes open at the new, reconnected Fermi surfaces when 
the hole and electron pockets fold down in the SDW state, despite the fact that the s± gap changes sign between 
the two pockets. Instead, the order parameter preserves its sign along the newly formed Fermi surfaces. The 
familiar experimental signatures of an s± symmetry are still preserved, although they appear in a mathematically 
different way. For a regular i case (s++) the nodes do appear in the SDW state. This distinction suggests a 
novel simple way to experimentally separate an s± state from a regular s in the pnictides. We argue that recently 
published thermal conductivity data in the coexisting state are consistent with the s± , but not the s++ state. 

PACS numbers: 74.20.Rp, 76.60.-k, 74.25.Nf, 71.55.-i 
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The superconducting pnictides continue to attract great 
interest over a year after the original discovery. Despite 
more than a thousand preprints and publications, the most 
basic questions about pairing symmetries and mechanisms 
remain controversial. Early on the s± pairing symmetry 
was proposed 12 , in which the superconducting gap function 
changes sign from the hole to the electron pockets, but is 
roughly constant on each surface. A possibility of an acci- 
dentally nodal s± state, or a d-wave state, depending upon pa- 
rameter values, was also proposed 2 , and investigated in many 
details recently 3 -. 

As of now, significant experimental evidence has been ac- 
cumulated in favor of the s± proposal, and substantial theoret- 
ical effort has been devoted to the study of various properties 
of such a state (see Refs. for reviews). So far, however, 
no one has addressed the possible modification of an s± state 
due to a static spin density wave (SDW) coexisting with su- 
perconductivity. At the same time the em erging consensus 
among experimentalists (see Refs. 617.8 9 BTolfTll) is that in 
most systems, most notably in both electron and hole doped 
122 materials, there is a range of coexistence of SDW and 
superconductivity, probably up to the optimal doping level 
(some, however, have argued for mesoscopic phase separation 
on the hole-doping side 11 ^). It was recently estimated that 
the magnetic moment at the Co concentration of 7% is 0. 1 [1b 
per Fe, corresponding, roughly, to an antiferromagnetic field 
of the order of 50-100 me ViS. 

The subject of an SDW coexisting with superconductivity 
has a long history, dating back to Bulaevskii et al in 1980 13 
and numerous work since then. It was show n 13 ' 14 that in a 
one-band BCS superconductor a spiral SDW induces a gap 
anisotropy that leads to gap nodes, while a collinear SDW 
still leads to a finite energy gap 14 . The case of the s± super- 
conductivity in Fe-based superconductors (FBS), on the first 
glance, seems quite simple: first, the SDW wave we are deal- 
ing with here is simply a collinear double-cell antiferromag- 
netic ( AF) order, so one need not be bothered by the difference 
between a spiral and collinear SDW; second, the doubling of 
the unit cell in real space leads to the folding down of the 
Brillouin zone in momentum space, which projects the elec- 



tron Fermi surfaces (FS) with the negative order parameter 
(A e < 0) onto the hole Fermi surface with the positive order 
parameter (A/, > 0). Whenever the two FSs intersect, an SDW 
gap opens up. It seems obvious that, when that happens, A 
on the newly formed FSs should change sign, that is, develop 
nodes. 

However, not everything that seems obvious is true. We 
will show below that, instead, a curious novel state is formed, 
which is fully gapped and, formally, has an order parameter 
(OP) of the same sign everywhere. This should not be con- 
fused though with the conventional BCS-like s state: when the 
SDW amplitude is vanishingly small, this state has the same 
observable properties as the original s± state, despite having 
a single-sign OP. This bizarre property, which, incidentally, 
is also relevant to the coexistence of d-wave superconductiv- 
ity and AF order in electron-doped cuprates 15 , can be traced 
down to two facts, well known but often not appreciated: (1) 
not only the overall sign of the OP in a superconductor, but 
also the relative sign of Au and A^i is not uniquely defined, 
but depends on the convention for the wave function phases 
and (2) as opposed to a nonmagnetic material, in an AF metal 
it is not possible to fix the phases of the wave functions in 
such a way that the wave functions for both spin projections 
are identical at any k-point. 

With these considerations in mind, let us now outline the 
derivation. We will follow the approach of Ref. [l4j, and 
for illustrative purposes will use a simple semimetallic model 
bandstructure with a hole band centered at the F point and an 
electron band around Q =(tc,0) and related points in the un- 
folded Brillouin zone which we will be using . We assume 
the Fermi energies to be, respectively, £/, and E e , and take an 
isotropic effective mass m for the hole band. To account for 
the fact that nesting is always imperfect (and if it were per- 
fect, the SDW would open a gap on the entire FS, thus pre- 
venting any coexisting superconductivity), we take m e to be 
anisotropic, with m x ^m y . This reflects the fact that the actual 
calculated and measured anisotropy of the electron pockets is 
larger than that of the hole ones. The Hamiltonian for this 
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system is given by 13 

M 

with an effective field /z °= cos(Qr), which interacts with the 
electron spin S and leads to the SDW. The index i = e,h refers 
to the hole or electron Fermi surface bands, with the disper- 
sions 

Z ke — h 2 (k x — %) 2 /2m x + h 2 k 2 /2m y — t e , 
£k.h = -tfk 2 /2m + Eh- 

For the purpose of this work - demonstrating the effect of AF 
- we will take £/, and e e equal. In real life, of course, they will 
depend on the relative location of Ep which will change with 
doping. We take m y <m< m y to ensure intersections between 
the hole and electron Fermi surfaces upon translating by the 
SDW vector. We take m x m y <~ m 2 so that the electron and hole 
densities of states (DOS) are comparable^ 

In the following we will work in the downfolded Brillouin 
zone corresponding to the antiferromagnetic unit cell. Let G 
be the matrix element (assumed to be k— independent) of the 
SDW potential mixing the hole and electron wave functions, 
cp/ 7 and (p,, . Then the dispersion in the SDW state is 

„ ± £k./, + £k,e ± v / ( e k./ ! -£k. e ) 2 +4G 2 

E k j (3) 

with the new wave functions 
\|/+ T = cose A .(p W! + sin9 A .(pi, e ; x\f^ = sin 9^.,/, - cosG^cp^^ 

(4) 

V||J = cos 8*9*,/, - sin9*9*, e ; \|/^ = sin9*9*,/, + cos9*9* jff 

(5) 

where tan 9k = G/(E^ — £k. e ). In Figure 1 we plot the sepa- 
rate hole and electron Fermi surfaces above the SDW ordering 
temperature (main panel), as well as the Fermi surface for two 
values of G. Of course, we are most interested in the limit in 
which the SDW-created reconstruction of the Fermi surface 
is relatively minor. In practical terms the SDW gap may be 
(albeit not necessary for retaining superconductivity!) smaller 
than the SC gap, but the concept is easier to illustrate for an 
SDW gap comparable to the superconducting one. Figure 2 
shows the behavior of Q k along the SDW-state Fermi surfaces. 
Note that while the energy spectrum (01 is double degenerate, 
the wave functions © are not equal to the the wave functions 
(0. If we now create a singlet anomalous average (the super- 
conducting OP), it will look like 

A + (k)= (y\fay± ki } = cos 2 9* (9^9^) (6) 

- sin 2 8* (q>k,e9k,e) + 2 sin 9* cos 8* (<Pk,e<?k,h) , 

and a similar expression for A~(k). At G — > 0, on the new 
"+" FS pockets 6* takes the value of or k/2 depending on 
whether the point on the FS originated from holes or from 
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FIG. 1: (Color online) Main figure: heavy line, the Fermi surface in 
the SDW state for G = EplYl\ dashed line, the Fermi surface above 
the SDW ordering temperature. The "+" and "— " Fermi surfaces are 
indicated. Inset: the Fermi surfaces for G = Ep / 120. 



electrons (and the reverse holds for the "— " pockets). Recall- 
ing the original s± assumption, A h = —A e with A'=(9,-9,-), 
(9/,9e) = we immediately observe that A + , as opposed to 
A h,e , never changes sign! 

However, the system remembers all too well that part of 
the new FS has come from electrons and part from holes, and 
that these used to have the OPs of the opposite sign. If we 
try to calculate any observable quantity (the OP per se is not 
observable), such as pair scattering from one part of the FS to 
another, we will have to take into account the fact that the new 
wave functions for the up spin are approximately equal to 
the old wave functions (e- or h-, depending on what part of 
the pocket we consider), but for the down spin © the same 
holds for the hole-type part of the FS, while for the electron- 
type parts the sign of the wave function is flipped. That is 
to say, any observable matrix elements include a product of 
the OPs and of the one-particle wave function; in the original, 
unfolded bands, we were able to choose signs of the wave 
functions to be spin-invariant, so that the scattering from the 
hole to the electron FS involved a sign change (s± ) . However, 
in the new, AF zone, no sign change is generated by the OPs, 
but the same sign change necessarily appears from the normal 
part of the matrix elements. 

We see that we cannot describe the system even with in- 
finitesimally weak SDW in terms of the same wave func- 
tions we used for the nonmagnetic parent system. Yet one 
can restore the conceptual continuity by describing the parent 
system differently. Let us select the wave function phases 
in the nonmagnetic system so that <$khl{ r ) — 9t,/i,t( r )' ^ ut 
9(fc,e,|( r ) = — 9(t,e,T( r )- This is a very inconvenient, but legit- 
imate gauge. In this gauge, the OP in the s± state will have 
the same sign on the both FSs, but any physical observable 
involving pair scattering will have to account for a sign flip 
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for electrons, but not holes, and thus the overall result will be 

unchanged — the same physical situation that we find (and in 

that case cannot avoid) in the SDW state. 

Let us move on to an arbitrary strength SDW and 

evaluate 17 the pairing matrix involved in a superconduct- 
or 

ing state below the SDW ordering temperature, A^, = 
-(V^lt/lit/^^^Jt/IV^), where U is the pairing interac- 
tion, with all relevant factors included and a, (} = ±1 are new 
band indices. We assume that only hole-electron matrix ele- 
ments {tyk,h\UWk' ,e) are nonzero, and the minus in front ac- 
counts for the fact that the pairing interaction is assumed to be 
generated by spin fluctuations. We furthermore assume that 
this matrix element does not depend on k, k! ^ 

After some trigonometric manipulations, we get the an- 
swer: 

A^, = ^(1 — apcos20£cos28 i / + apsin28 i .sin20 ; (./) 

The factor 1/2 was selected so that (as we will show later) the 
effective coupling constant in the G = limit will be equal to 
Vy/N e Nh. Now we can write the BCS equation at T = T c as 
(with co being the BCS cutoff energy in temperature units) 



ln(1.13to/T c )A« = ^£8(4)A, 



2 *— 

£ k'R 



(7) 



V 



-acos 2Q k £ P5(£[, ) cos 2Q k , A k , 



/ ( 'p 



-asin26^p8(4)sin29^Af 



(8) 



We seek the solution of this equation in the following form 



A^/Aq =c + aacos2e i + a/7sin2e A: , (9) 



and 



ln( 1 . 13<d/T c )c = - (aN c + bN s + cN) 

ln(1.13co/r c )fl = --(aN a +bN s + cN c ) (10) 

ln(1.13co/T c )Z> = -(aNa + bNa + cN,} 

where we introduced the weighted DOSs: N = £8 (£■"), 
N c = Ea8(£?)cos29i, N s = ZaS(E?)sm2Q k , N c2 = 
l5(£f)cos 2 2ei, N s2 =N-N c2 = I8(£«)sin 2 29,t, N cs = 
£8(£'")cos20 / ( : sin28< : . All summations are overfc and a, and 
8 represents the Dirac delta function. 
The maximal eigenvalue of the matrix 



V 



N N c N s 
-Nc ~Nc2 ~N CS 
N N cs N s2 



(11) 



as usual, defines the BCS transition temperature T c , and the 
corresponding eigenvector gives the distribution of the order 
parameter over the Fermi surfaces. 



In the limit G -> 0, cos29 i: -> ±1, N c -> N e -N h ,N c2 -> N, 
and all other DOSs vanish. Here N e ^is the density of states 
on the electron (hole) Fermi surface without th e SDW. The ef- 
fective coupling constant X e ff = V ^N 2 — /2 = Vy/N e Nh, 
which is the well-known result for the s± pairing. The ratio of 
the two gaps at G — > is ^Nh/N e as it should be in a weak- 
coupling s± superconductor^ 

At a finite G, the gap does become angle-dependent, via 0^, 
but it is easy to proved that the eigenvector for maximal A, 
in © gives c > a ^S> b, so that the gap has the same sign ev- 
erywhere, in agreement with general discussion after Eq. (O. 
The gap is depicted in Figure 2, where we note that the gap 
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FIG. 2: (Color online) Depicted is the variation of A(k) (solid lines) 
and Q k (dashed lines) along the Fermi surfaces for the indicated val- 
ues of G. Note that the angle 9^ will approach a step function, with 
regions where Q k = or 9^ = 7l/2, as G — > 0. Here "length" refers to 
the k-space arc-length around the Fermi surface; note that the actual 
arc-length depends on G. To facilitate comparison we have scaled the 
arclengths for these two cases lengths to be equal; actual G = Ef/6 
arclengths are ~ 20 percent smaller than for G = Ep f 12. 

changes rapidly within a finite area where the SDW gap opens; 
the length of this region scales with G. 

We have also found that the effective coupling constant is 
reduced slightly at a finite G compared to its value at G = 0, 
i.e., superconducting T c drops in the presence of an SDW 
order. At small G, we found 8X e ff ~ — (VN s2 /%)(\/N — 
VN~h) 2 /VNM and 20 N s2 < 



Note that the de- 



\G\/\N e -N h \ 
pendence on G is non-analytic. 

The calculations above can be easily generalized to three 
other cases: an s± with a charge density wave (CDW), a regu- 
lar s with a SDW, and a regular s with a CDW. The last case is 
mathematically equivalent to the one considered above. The 
first two cases are equivalent to each other and constitute a dif- 
ferent set. Indeed, in that case the V in Eqs. 171 [Tol [TTI should 
have the opposite sign, and the solution will always have eight 
node lines at the tips of the banana-shaped FSs of Fig. 1 . The 
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eigenvector for the largest eigenvalue will have, for N e ~ N e , 
the largest weight on a, not c, and the OP will have one sign 
for the parts of the FS that originated from the electrons and 
the opposite sign for those originating from the holes. 

To summarize, we have shown that, surprisingly, the SDW 
observed in underdoped pnictide compounds, does not have 
any considerable destructive effect on the s± superconductiv- 
ity, besides the obvious competition between the two instabili- 
ties for the density of states at the Fermi level. As opposed to a 
hypothetical CDW, which would have created nodes on the FS 
and additionally weakened superconductivity, an SDW wave 
retains the gapped nature of the s± superconductivity. Never- 
theless, this constant-sign state has the same observable phys- 
ical properties as the sign-changing s± state wit hout an SDW. 
In particular, the penetration depth X(T) °< 1 /y/ps(T) is still 
exponential at the lowest T and crosses over to roughly T 2 
at higher T-. The slope of SK(T) °c 8p s (T)/(fi s {T = 0)) 3 / 2 
however increases with increasing SDW order simply because 
p s (T = 0) decreases together with the area of the Fermi sur- 
face. This slope increase upon approaching the SDW state 
was observed in Ref. [22l If, however, the pairing state were 



a conventional s-wave, an SDW order would give rise to gap 
nodes, and 8A- would be linear in T at the smallest T (we have 
explicitly verified that the SDW coherence factors do not af- 
fect this behavior). Similarly, the reduced thermal conduc- 
tivity k/T should still vanish at T — ► if the gap is s± but 
should become finite below the onset of the SDW order if the 
gap is a conventional s-wave, and the more we go into the 
underdoped regime, the larger the residual k/T should be be- 
cause its value does not depend on impurity concentration, 
but only on the inverse of the slope of near the nodes. In- 
deed, a recent study of the low-T thermal conductivity in the 
underdoped BaFe2-. v Co x As2, where microscopic coexistence 
of SDW and superconductivity has been well documented, 
shows 23 that k/T vanishes at T — > indicating absence of any 
gap nodes. This is consistent, according to our results, with an 
s± state. 
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